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Recently, Hofava proposed a power counting renormalizable theory for (3+l)-dimensional quan- 
tum gravity, which reduces to Einstein gravity with a non-vanishing cosmological constant in IR, 
but possesses improved UV behaviors. In this work, we analyze the stability of the Einstein static 
universe by considering linear homogeneous perturbations in the context of an IR modification of 
on : Hofava gravity, which implies a 'soft' breaking of the 'detailed balance' condition. The stability 

regions of the Einstein static universe is parameterized by the linear equation of state parameter 
w — p/p and the parameters appearing in the Hofava theory, and it is shown that a large class of 
stable solutions exists in the respective parameter space. 
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Recently, motivated by the Lifshitz model in condensed matter physics, Hofava proposed a power counting renor- 
i 1 malizable theory for (3+l)-dimensional quantum gravity [l|, 0]. This theory, denoted as Hofava-Lifshitz gravity, is 
believed to be the potential ultraviolet (UV) completion of general relativity (GR). In the infrared (IR) limit (setting 
the parameter A = 1 in the action), it recovers GR. Hofava-Lifshitz gravity admits a Lifshitz scale- invariance in time 
and space, exhibiting a broken Lorentz symmetry at short scales, while at large distances higher derivative terms 
do not contribute, and the theory reduces to standard GR. Since then various properties and characteristics of the 
Hofava gravities have been extensively analyzed, ranging from formal developments Q, cosmology Q, dark energy 
^C) ■ @> @ an d dark matter Q , spherically symmetric solutions [H, [^] , and its viability with observational constraints flQ ] 
were also explored. Although a generic vacuum of the theory is the anti-de Sitter one, particular limits of the theory 
allow for the Minkowski vacuum. In this limit post-Newtonian coefficients coincide with those of pure GR. Thus, the 
deviations from conventional GR can be tested only beyond the post-Newtonian corrections, that is for a system with 
| strong gravity at astrophysical scales. 

In this work, we consider the stability of the Einstein static universe in Hofava-Lifshitz gravity, with a 'soft' 
violation of the detailed balance condition (Recently, the stability of the Einstein static universe in Hofava-Lifshitz 
gravity satisfying the detailed balance condition was analyzed (HI). The presence of the respective term in the action 
which represents a 'soft' violation of the 'detailed balance' condition modifies the IR behavior. Note that this IR 
modification term, with an arbitrary cosmological constant, represent the analogs of the standard Schwarzschild- 
(A)dS solutions, which were absent in the original Hofava model. The analysis of the static Einstein Universe is 
motivated by the possibility that the universe might have started out in an asymptotically Einstein static state, in 
the inflationary universe context [l|. On the other hand, the Einstein cosmos has always been of great interest in 
various gravitational theories. 

In GR for instance, generalizations with non-constant pressure have been analyzed in jl3| . In brane world models, 
the Einstein static universe was investigated in [T3 |. while its generalization within Einstein-Cartan theory can be 
found in [TEl ]. and in loop quantum cosmology, we refer the reader to [jq |. In the context of f(R) modified theories of 
gravity, the stability of the Einstein static universe was also analyzed by considering homogeneous perturbations [l7| . 
By considering specific forms of f(R), the stability regions of the solutions were parameterized by a linear equation 
of state parameter w = p/p. Contrary to classical GR, it was found that in f(R) gravity a stable Einstein cosmos 
with a positive cosmological constant does indeed exist. Thus, in principle, modifications in f(R) gravity stabilize 
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solutions which are unstable in GR. Furthermore, in [18] it was found that only one class of f(R) theories admits 
an Einstein static model, and that this class is neutrally stable with respect to vector and tensor perturbations for 
all equations of state on all scales. These results are apparently contradictory with those of Ref. [lj]- However, in a 
recent work, homogeneous and inhomogeneous scalar perturbations in the Einstein static solutions were analyzed [l9| , 
consequently reconciling both of the above works. In the context of modified theories of gravity, the stability of the 
Einstein static universe in /(G) Gauss-Bonnet modified gravity was also analyzed [20l |. In particular, by considering 
a generic form of /(G), the stability regions of the Einstein static universe were parameterized by the linear equation 
of state and the second derivative /"(G) of the Gauss-Bonnet term. It was shown that stable modes for all equation 
of state parameters w exist, if the parameters of the theory are chosen appropriately. Thus, the results show that 
perturbation theory of modified theories of gravity present a richer stability /instability structure than in GR. 

Thus, it is the purpose of the present paper to consider the stability of the Einstein static universe by considering 
linear homogeneous perturbations in Hofava-Lifshitz gravity. Indeed, this analysis is particularly important as the 
higher derivative terms in the action contributes with a 1/a 4 term in the modified Friedman equations. This contri- 
bution becomes dominant for small a, and as mentioned above motivates this analysis, due to the possibility that the 
universe might have started out in an asymptotically Einstein static state, in the inflationary universe context [l2l ]. 
On the other hand, the cosmological solutions of GR are recovered at large scales. It is shown that a large class of 
Einstein static universes exist that are stable with respect to linear homogeneous perturbations. 

This paper is outlined in the following manner: In Sec. |TIJ we briefly review the action and field equations of Hof ava 
gravity, and the respective modified Friedman equations. In Sec. IIII1 we consider linear homogeneous perturbations 
in the context of the Einstein static Universe in Hofava gravity, and analyze the respective stability regions. In Sec. 
HVl we conlcude. 

II. HORAVA GRAVITY AND FIELD EQUATIONS 
A. Action 

Using the ADM formalism, the four-dimensional metric is parameterized by the following 

ds 2 = -N 2 c 2 dt 2 + gij (dx l + N* dt) (dx j + N 3 dt) , (1) 

where N is the lapse function, N' 1 is the shift vector, and is the 3-dimensional spatial metric. 
In this context, the Einstein-Hilbert action is given by 

S = ^ J d 4 x^N (r^K* -K 2 + - 2A) , (2) 

where G is Newton's constant, is the three-dimensional curvature scalar for g^, and K%j is the extrinsic curvature 
defined as 



where the overdot denotes a derivative with respect to t, and V, is the covariant derivative with respect to the spatial 
metric g^. 

Consider the IR-modified Hofava action given by 



2 

^2 



Kij K- - XK 2 ) - ^ Cij C« + ^ k Rfv 3 R^ k 



S = J dtd 3 t ^/g~N 



_ R (3) 



(4) 



8(3A- 1) 

where K, A, v, /1, vo and A^y are constant parameters. C lJ is the Cotton tensor, defined as 

cij = £lMVfc f R (3)j i _ l R (3) 5 A . (5) 



Note that the last term in Eq. (01 represents a 'soft' violation of the 'detailed balance' condition, which modifies the 
IR behavior. This IR modification term, /j, 4 R^ 3 \ generalizes the original Hofava model (we have used the notation of 
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Ref. 0). Note that now these solutions with an arbitrary cosmological constant represent the analogs of the standard 
Schwarzschild-(A)dS solutions, which were absent in the original Hofava model [9(. 

The fundamental constants of the speed of light c, Newton's constant G, and the cosmological constant A are 
defined as 



K 2 p 2 \A w \ 
8(3A- l) 2 



G = 



K 2 C 2 



16tt(3A - 1) 



A = -A w c 2 . 



(6) 



B. Modified Friedman equations 



Consider the homogeneous and isotropic cosmological solution given by the following metric 

dr 2 



ds 2 



-dt 2 + a 2 (t) 



1 - kr 2 



+ r 2 (d6 2 + wo* 6 dtp) 



(7) 



where k = +1, 0, —1 corresponds to a closed, flat, and open universe, respectively. 

We assume that the matter contribution takes the form of a perfect fluid, with p and p the energy density and the 
pressure, respectively, so that the modified Friedman equations in Hofava gravity take the following form 
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6(3A- 1) 
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8(3A- 1) V a 
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(8) 
(9) 



where e = ±1 . 

The analytic continuation p 2 — > — p? for the dS case, i.e., A\y > 0, is considered @, and the upper (lower) sign 
denotes the AdS (dS) case. It is interesting to note that the higher derivative term appearing in the action Q 
contributes with a 1/a 4 term, and only exists for k ^ 0. This term dominates for low values of a, and the general 
relativistic cosmological solutions are recovered for large scales. 



III. THE EINSTEIN STATIC UNIVERSE IN HORAVA GRAVITY AND PERTURBATIONS 

A. Field equations 

For the Einstein static universe, a = fto = const and k = 1, the Ricci scalar becomes R = 6/a 2 ,, (note that for 
A = 1, GR is obtained in the IR limit). Furthermore, we consider a linear equation of state, p = wp, so that the field 
equations in this case are expressed in the following manner 



A) 



2 ? 
€K (J, 



2 + 3a§(l + w)vo ± s/4 - 6a 2 , (1 + w)(l + 3wjza 



6(l + w) 2 (3A-l) 



(10) 



2 (1 + 3w) T V4 - 6a 2 (l + w)(l + 3w)ro 
3(1 + w) 

where po and po are the unperturbed energy density and isotropic pressure, respectively. Note that it is useful to 
introduce the dimensionless parameters fl := a^zu and A = AwOq- These relationships are useful to be written in 
this form, as one has a first glance at the existence issue that po and A should be real, which imposes the following 
condition 

2-3(l + w)(l + 3w)f2>0. (12) 

The above inequality needs to be analyzed for the three cases w < —1, —1 < w < —1/3 and w > —1/3 which places 
restrictions on the allowed values of Q. However, as we are primarily interested in a physically reasonable Einstein 
static universe, we will furthermore require positivity of the energy density, po > 0. To further simplify the subsequent 
analysis, we will also assume A > 1/3. 
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These useful conditions imply the following existence conditions for the upper sign of Eq. ([10 

2 1 

e = +1 , w < -1, <C — 



3(l + u))(l + 3to) ' 

+1, -l<u><-2/3, n>| 7 — - — , 

3 (1 + u>)(l + 3u>) 

+1, - 2/3 < w < -1/3, 0>-2, 

2 1 

+ 1, -l/3<w, -2<fi<-- 



3 (1 + w)(l + 3w) 

-I, -2/3<™<-l/3, H J_ <^<-2, 

3 (1 + w)(l + 3w) 

-1, -1/3 <w, n<-2, (13) 



which are depicted in the left plots of Figs. [T] and [21 respectively. 
For the lower sign, of Eq. (flO)) , we find: 

e = +1, w < -1, < -2, 

2 1 

e = +1 , w < -1 , < 17 < - 



3 (l+w)(l+3w) ' 

-1, -l<«;<-2/3, f — J <0<-2, 

3 (1 + w){l + 6w) 

-1, -1<w<-1/3, < O , 

2 1 

-1, -l/3<tu, < S7 < 



3(l + w)(l + 3w) ' 
e = -l, it; < -2/3, -2<O<0, 

^"^ -V3<«<-l/3, § (1 + tl>) ( 1 + 3a>) <"<0, 

e=-l, -l/3<w, il <0. (14) 

which are depicted in the right plots of Figs. [T] and [21 respectively. 

The above inequalities with the upper and lower signs, and with e = +1, are represented graphically in Fig. [U 
and with e = — 1 in Fig. [21 respectively. Note that there are many configurations which allow for an Einstein static 
universe. This is in contrast with other modifications of GR, such as f(R) modified gravity, where there exists a 
unique background Einstein static universe for every choice of f(R)- 

B. Linear homogeneous perturbations 

In what follows, we analyze the stability against linear homogeneous perturbations around the Einstein static 
universe given in Eqs. Thus, we introduce perturbations in the energy density and the metric scale factor 

which only depend on time 

p(t) = p + Spi(t), a{t) = a + 6a 1 (t). (15) 

Now we consider adiabatic perturbations which also satisfy a linear equation of state, Sp(t) — w5p(t) and linearize 
the perturbed field equations. Firstly, we consider Eq. (jSJ) which upon subtracting the background field equation 
yields 



3en 2 p? 
2ajj(3A- 1)' 



frl= o,5,ox ^ (A- 0-1) fax. (16) 



Next, we perturb the evolution equation @ and eliminate the perturbed energy density by virtue of the latter 
equation. The resulting second order differential equation for 5a\ is given by 



4, ,2 



K(t) + 8a 4 (3A - IV [(1 + A " Q) ~ 3W{1 ~ A + * ai (t) = ° • 



8a 4 (3A - 1) 
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Existence regions Existence regions 




FIG. 1: Regions of existence in the (w, fi) parameter space for the specific case of e = +1. The left panel shows the existence 
regions for the upper sign of Eq. (110[) , while the right panel is for the lower sign. The different shades of gray correspond to 
the different inequalitites. 



Existence regions Existence regions 




FIG. 2: Regions of existence in the (w, SI) parameter space for the specific case of e = — 1. The left panel shows the existence 
regions for the upper sign of Eq. (|10p . while the right panel is for the lower sign. The different shades of gray correspond to 
the different inequalitites. 



As only the sign of the prefactor of the second term is relevant, we can rescale 8a\(t) appropriately and consider 

5a"(t) + e[(l+A-Q) + 3w(-l + A-Q)]5ai(t) = 0. (18) 
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Before solving this equation it should be noted that A is determined by the background, see Eq. (fTTj) and hence this 
quantity should also be substituted. We find 



(1 - 3ft) + 3to(l - ft) + V 4 - 6(1 + w)(l+ 3w)Q , . 

A-ft = - i ' v r — . (19) 

3(1 + to) V ' 

Using the standard ansatz a\{t) — Aexp(iWt), where A and W are constants, we find that the above differential 
equation (|18p provides the following solutions 



W = ±y/ey/(l + A - ft) + 3w(-l + A - ft) . (20) 



C. Stability regions 

To analyze the stability regions, first consider the cases e = +1 (e = — 1). Now, the stability conditions impose that 
the factor within the square root of Eq. (|20[) is positive (negative), which is translated by the following inequalities: 

e = +l (l + A-ft)+3iu(-l + A-ft) >0, (21) 
e = -l (l + A-ft)+3iu(-l + A-ft) <0. (22) 

Since A is determined by the background equations, see (|19p . the equations become more complicated. For the 
upper sign we have: 

2 1 

e = +i , - 1 < io < -i/3, ft > 



3(l + w)(l + 3w) 



e = +l, -1/3<iu<1/3, ft<^- ' 



3(l + w)(l + 3w) 



_ l-3to 1 -l + 3w 

e=+l, l/3<w, 0< T^ W-n 

1 + 3w V3 V 1 + w 



l-3u; 1 /-l + 3w rt 2 
-1, 1/3 <w, — + -=J— <ft<-- 



l + 3w 73 V l + w 3(l + w)(l + 3w) ' 



2 1 

e = -1 , to < -1 , ft < 77- 



3(l + w)(l + 3w) 



l-3io 1 /— l + 3*p _ l-3w 1 /-l + 3w 
1 + 3w V3v l + w 1 + 3w V3v l + w 

which are depicted in the left plots of Figs. OandlH respectively. 
Finally, for the lower sign we obtain: 



l-3w 1 -l + 3w _ l-3w 1 -l + 3w 

e = +1 ' w< - 1; T^-vlVTT^r < ^ < TT3w + v!VTT^r 

2 1 

e = +l, -l/3<to, ft < - 

3 (1 + w){l + 3wj 



„ l-3w 1 -l + 3w 

w<-l, ft< = 4/ 

1 + 3w V3 V l + w 



!- 3u; 1 /-l + 3w _ 2 1 

e = -l, w<-l, — =\ < ft < —. -, (24) 

l + 3w V3V l + w 3 (1 +w)(l + 3w) ' v 7 

which are depicted in the right plots of Figs. [3] and |3J respectively. 

By combining the existence conditions with the stability conditions, we can identify a large class of Einstein static 
universes which are stable with respect to homogeneous perturbations. It should also be noted that for the upper sign 
a Einstein static universe of phantom matter (w < — 1) exists. Superimposing the inequality plots Figs.[lHl]we can 
picture the complete parameter space for which the Einstein static universe in IR modified Hofava gravity exits and 
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FIG. 3: Regions of stability in the (to, Q) parameter space (left panel, upper sign; right panel, lower sign), for the specific case 
of e = +1. See the text for details. 
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FIG. 4: Regions of stability in the (w,Q) parameter space (left panel, upper sign; right panel, lower sign), for the specific case 
of e = — 1. See the text for details. 



is stable with respect to linear homogeneous perturbations. This is depicted in Fig. [5] for the specific case of e = +1. 
For the specific case of e = — 1, we verify the nonexistence of the Einstein static universe existence/stability regions. 

For every parameter choice we can compute the actual value of A by using Eq. (fTTjl . Note, however, that it is very 
involved to represent the values A can attain in general because this would require to plot a surface in the (to, f2, A) 
parameter space which satisfies various inequalities simultaneously. 




IV. SUMMARY AND DISCUSSION 



The Einstein static universe has recently been revived as the asymptotic origin of an emergent universe, namely, as 
an inflationary cosmology without a singularity [13]. The role of positive curvature, negligible at late times, is crucial 
in the early universe, as it allows these cosmologies to inflate and later reheat to a hot big-bang epoch. An attractive 
feature of these cosmological models is the absence of a singularity, of an 'initial time', of the horizon problem, and 
the quantum regime can even be avoided. Furthermore, the Einstein static universe was found to be neutrally stable 
against inhomogeneous linear vector and tensor perturbations, and against scalar density perturbations provided that 
the speed of sound satisfies c|? > 1/5 [23(. Further issues related to the stability of the Einstein static universe may 
be found in Ref. H. 

In this work we have analyzed linear homogeneous perturbations around the Einstein static universe in the context 
of IR modified Hofava gravity. In particular, perturbations in the energy density and the metric scale factor were 
introduced, a linear equation of state, p(t) — wp(t), was considered, and finally the linearized perturbed field equations 
and the dynamics of the solutions were analyzed. It was shown that stable modes for all equation of state parameters 
w exist and, in particular, the complete parameter space for which the Einstein static universe in IR modified 
Hofava gravity exits and is stable with respect to linear homogeneous perturbations was presented. Thus, as in 
Refs. pj, H, our results show that perturbation theory of modified theories of gravity present a richer 
stability/instability structure that in general relativity. Finally, it is of interest to extend our results to inhomogeneous 
perturbations in the spirit of Ref. [l9| , and to include the canonical scalar field case. Work along these lines is presently 
underway. 
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